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I briefly review lattice QCD, the valence (quenched) approximation, and the ap- 
plication of both to the determination of the mass and two-body decay couplings 
of the lightest scalar glueball. Results in agreement with the observed properties 
of /j(1710) strongly suggest this resonance is largely a scalar glueball. 

1 Introduction 

The existence of chromoelectric field is the key hypothesis QCD adds to the 
quark model. The observation in experiment of a glueball would be a direct 
confirmation of the existence of chromoelectric field. So it would be nice to find 
one. Identifying a glueball in experiment is tough, however. The kinematics of 
QCD provides no clear footprint for glueballs comparable, say, to the footprint 
the Weinberg-Salam model gives for the W's and the Z. The lightest states, 
which would be easiest to detect, are not expected to differ drastically in mass 
or decay properties from mesons containing quarks and antiquarks. On the 
other hand, J PC combinations which are impossible for quark-antiquark states 
but which can occur for heavier glueball excitations can also be realized with 
multiple quark-antiquark pairs. A reliable calculation of the consequences of 
QCD's dynamics for the masses and decay properties of glueballs appears to be 
a necessary element in the identification of glueballs in experiment. At present 
such calculations can be done only numerically, using the lattice formulation 
of QCD in combination with the valence approximation. 

As of three years ago, two independent calculations had been completed 
of the infinite- volume, continuum limit of the mass of the lightest glueball, 
which turns out to be_a scalar. Hong Chen, Jim Sexton, Alessandro Vac- 
carmo and I obtained tH a value of 1740(71) MeV using ensembles of 25000 
to 30000 gauge configurations on each of several different lattices. An earlier 
valence approximation calculation by the JJKQCD-Wuppertal collaboration^, 
when extrapolated to the continuum limitLl, yields 1625(94) MeV for the light- 
est scalar glueball mass. This calculation used several different lattices with 
ensembles of between 1000 and 3000 configurations each. If the two mass 
evaluations are combined, taking into account the correlations between their 
statistical uncertainties arising from a common procedure for converting lattice 
quantities into physical units, the result is 1707(64) MeV for the scalar glueball 
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mass. Both the mass prediction with larger statistical weight and the combined 
mass prediction agree with the mass of /j(1710) and are inconsistent with all 
but /o(1500) among the established flavor singlet resonances which could be 
scalars. For /o(1500) the disagreement is still by more than three standard 
deviations. In addition, observed mass values for other scalar quark-antiquark 
states suggest a mass near 1500 MeV for the lightest ss scalar. 

These calculations by themselves, however, do not make a strong case for 
taking /j(1710) as a scalar glueball. The key question which the mass results 
do not answer is whether the lightest scalar glueball has a decay width small 
enough for this particle actually to be identified in experiment. If the scalar 
glueball had a width of a GeV or more, the prospect of ever finding one seems 
remote. But a scalar glueball with a width of a few hundred MeV or less and 
mass in the neighborhood of 1700 MeV would be hard to miss and should 
already have been seen in experiment. A further question in the identification 
of /,/(1710) as a glueball is raised by the argument that since glueballs are 
flavor singlets they should have the same decay couplings to ttqtto, to K$Ks, 
and to r/r/. This equality is violated by /j(1710) decay couplings. 

To address these questions, Jim Sexton, Alessandro Vaccarino and 10 have 
calculated the decay coupling constants of the lightest scalar glueball to pairs of 
pseudoscalar mesons. The calculation was done in the valence approximation 
on a lattice 16 3 x 24 with /3, defined as G/g 2 , set at 5.70, corresponding to lattice 
spacing, determined from the p mass, of about 0.15 fm and lattice period of 
about 2.3 fm. For the total width of the scalar glueball to pairs of pseudoscalar 
quark-antiquark states we obtained 108 ±29 MeV. The combined correction for 
the errors in our prediction arising from the valence approximation, from finite 
lattice spacing and from finite lattice volume we believe would be less than 
50%. Based on our value for two-body decays, any reasonable guess for the 
partial width for multibody decay modes leads to a total width small enough 
for the scalar glueball to be hard to miss in experiment. In fact, the observedB 
width of /j(1710) into pairs of pseudoscalars is 99 ± 15 MeV, consistent with 
our result. We obtain also a violation of the expected equality of glueball decay 
rates to 7ro7ro, to K$Ks, and to r\r\ in agreement with the observed results for 
/./(1710). 

Although so far I have simplified the story by supposing that physical 
resonances are either entirely glueballs or entirely quark-antiquark, another 
possibility is a state which is a linear combination of a glueball and a quark- 
antiquark state. In the valence approximation, however, glueballs contain no 
admixture of configurations with valence quarks or antiquarks. Thus we con- 
sider the agreement between the mass and decay couplings found in the valence 
approximation and the observed mass and decay couplings of /j(1710) to be 
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strong evidence that this state is largely a scalar glueball with at most some rel- 
atively smaller amplitude for configurations including valence quark-antiquark 
pairs. 

The glueball calculations presented here were carried out on the GF11 
parallel computer B at IBM Research and took approximately 30 months to 
complete at a sustained computation rate of between 6 and 7 Gflops. 

In the remainder of this talk, I will briefly review lattice QCD, the valence 
approximation, the glueball mass calculation, the decay calculation, and close 
with a more detailed comparison of the lattice predictions with experiment. 



2 Lattice QCD 

Lattice QCD, as its name implies, approximates continuous space-time (at 
negative imaginary values of time) by a discrete lattice including only a finite 
number of points. Predictions for the real world of continuous, infinite-volume 
space-time are supposed to be recovered from lattice QCD by taking the limit of 
lattice predictions as lattice spacing goes to zero and lattice volume to infinity. 
It is convenient, though not necessary, to take the lattice to be hypercubic 
with periodic boundary conditions. Living at each site are lattice versions of 
the chromoelectric potential, A>(x), given by real numbers, and the quark 
and antiquark fields, ^S> s d{x) and ^ sc (x), respectively, given by Grassmann 
variables. Here x is a lattice site, j, c and d are color indices for the adjoint, 
fundamental and conjugate representations, respectively, fi is a lattice direction 
and s is a spin index. 

Vacuum expectation values of time ordered products of fields are then 
defined by a path integral. A typical vacuum expectation value, which I offer 
in place of the slightly more complicated general rule, is 

<*.).(.)> - r '/n*«/n****w***+*>. 

Z = JUdfiA jY[d^d^exp{S A +S^), (1) 



where J Yl d^A can be thought of, but is not quite, a product over all x, /i 
and j of the integral dAj^ (x) and / Y\ d^d^S is a product of Grassmann 
integration on each ^ S( i(x) and ty sc (x). For the continuum quantities and 
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integrals which appear in Eqs. ([j]), I actually intend lattice approximations 
but use the continuum expressions as more recognizable alias's. 

With respect to the quark and antiquark fields, Eqs. ([!]) are a Grassmann 
version of gaussian integrals. The integrals can be done analytically and give 

< ¥(x)\P(x) > 



Z 

where the trace is with respect to spin and color indices. 

The lattice QCD industry does integrals like those in Eq. @ by Monte 
Carlo. With a large enough ensemble of random A fields, [A^(x)]i, . . . [A^(x)]jv, 
generated by computer according to the differential probability 

dv = Z^ 1 dfiA det(j> — A 1 + im) exp(<S^). (3) 

the vacuum expectation value of Eqs. (|J) becomes 

<*(x)*(x) >=—J2 tr [(t-4k+im)- 1 {x,x)} (4) 
fe 

3 The Valence Approximation 

The problem in generating, with present computer hardware, ensembles of A 
according to dv of Eq. (||) , and the origin of the valence approximation, is the 
factor det($ — A 1 + im). For the moment forget about this factor. A simple 
algorithm to generate an ensemble of A begins with some arbitrarily chosen 
field, for example all A components set to 0, and then successively walks across 
all sites x and directions /i on the lattice modifying the vector of A^ (x) , for all 
j, at each. The collection of A resulting from each complete sweep updating 
all x and /i gives a possible Monte Carlo ensemble. The update at each x 
and /x proceeds in two steps. First a random new trial vector of Aj^(x)' , for 
all j, is generated according to a rule the details of which I will skip. Then 
the vector of A^(x)' is either installed in place of the old A^(x) or thrown 
away with a probability depending on the change AS a which this replacement 
would cause in Sa ■ Now, the action Sa approximates the derivatives in (x) 
with nearest-neighbor differences. So computing the change ASa arising from 
the trial replacement of A^{x) by A^(x)' leads to work involving only A*(y) 
for y equal to x or nearest-neighbors of x. Only a fixed number of arithmetic 



Z 1 J Y]_dfiA tr[(j/> - 4 + im)~ 1 (x,x)} x 
det(j/> — A" + im) cxp(Sa), (2) 
J IT d e t{i > ~ 4- + * TO ) ex P(^A), 
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operations are required independent of the size of the lattice. The total work 
to update (x) for all x and [i to produce a new member of the Monte Carlo 
ensemble is therefore proportional to V, the number of sites in the lattices. 

Now put det{j>— Jft+im) back in dv in Eq. (H). The algorithm for generating 
a Monte Carlo ensemble is nearly unchanged, except that in place of the the 
change ASa, we need the change ASa + A \og{det(]f) — A 1 + im)]. Since det{^> — 
A 1 + im) couples together all A^{x) on the entire lattice, the work to find 
ASa + A\og[det{]/) — jJl + im)] in the update for a single x and fi is already 
of order V. The work to update all sites becomes of order V 2 . But recall, 
for a moment, that predictions for the real world are found from lattice QCD 
predictions by taking limits of zero lattice spacing and infinite lattice volume. 
Even a rough approximation to these limits, at least for the simplest lattice 
implementations of Sa and S\^ in Eqs. (Q), requires lattice dimensions of at least 
10 x 10 x 10 x 10 thus V of 10 4 . The cost of including det(j/> - 4 + im) therefore 
becomes a factor of 10 4 or more. Actually, it turns out that somewhat fancier 
algorithms than what I just described handle det(^> — A 1 + im) at lower cost, 
perhaps a factor of 100 to 1000 for the largest lattices. Unfortunately, present 
computer power is just barely sufficient to handle lattices large enough to give 
infinite- volume, continuum limit predictions in the absence of det(fi — 4 + im) 
in dv. The increase by a factor of 100 to 1000 in work required to include 
det(j} — A + im) in dv is fatal. 

n 

What to do? A possible answer El is suggested by viewing det(j> — 4 + * m ) 
through the eyes of weak coupling perturbation theory. In weak coupling 
perturbation expansions, det(fi — jjL + im) gives rise to closed quark loops inter- 
rupting gluon lines inside diagrams. So det{$ — Jji+im) is in charge of the QCD 
analogue of the particle-hole polarization process that occurs when an electro- 
magnetic field propagates through a solid. In the case of electrodynamics, we 
know that for sufficiently weak, sufficiently low momentum electromagnetic 
fields the effect of particle-hole polarization is accurately approximated by 
simply replacing all charges and fields by screened values obtained by dividing 
each by the solid's dielectric constant. A plausible hypothesis is that the same 
holds for QCD. Namely, for processes involving sufficiently weak, sufficiently 
low momentum chromoelectric fields, we can omit det{j> — jjL + im) from dv 
and just replace the chromoelectric charge g in Sa of Eqs. (|l|) with a screened 
charge g/rj, for some QCD analogue dielectric constant 77. This replacement is 
the valence approximation, also called the quenched approximation. 

In part as a test of the method, a calculation! has been done of the infinite- 
volume, continuum limit of the valence approximation to light hadron masses. 
A total of 11 masses were calculated, three of which were used to determine 
input parameters. The ir mass was used to set the up and down quark masses, 
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Figure 1: Continuum limit valence approximation hadron mass predictions in finite volume, 
triangles, and extrapolated to infinite volume, circles. 

taken to be equal. The K mass was used to set the strange quark mass. The 
p mass was used to determine the gauge coupling constant. QCD's supply of 
free parameters is then exhausted. The remaining eight masses are predictions. 
The results in units of the p mass are shown in Figure [l]. The triangles give zero- 
lattice-spacing results in a box with period of about 2.3 fm, which turns out 
to be nearly infinite volume for light hadron masses. The circles show infinite- 
volume predictions obtained by applying an additional correction to the 2.3 
fm results. Out of eight numbers, the biggest disagreement between prediction 
and experiment is 6%. The statistical uncertainties in these numbers range up 
to 8%, however. For eight numbers with these uncertainties, one disagreement 
with experiment by 6% is expected. Thus 6% should be treated as the one 
sigma upper bound on the error in the valence approximation. 

So taking the valence approximation as accurate to within 6% for light 
hadron masses, we applied it to predicting glueball properties. 



4 Scalar Glueball Mass 



To determine the scalar glueball mass El we evaluated the vacuum expectation 
value 



C(t) = <g(t)g(0)>, 

g(t) = h(t)-<h(t)>, (5) 
h{t) = J <PxFi b (x,t)Fi b (x,t), 
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Figure 2: Scalar glueball effective mass on a lattice 30 X 32 2 X 40 with lattice spacing of 

0.050 fm. 

where a and b are summed only over space directions, and as before, I actually 
intend lattice approximations for F 3 ab (x,t) and J d 3 x but use the continuum 
expressions as more recognizable alias's. By inserting a complete set of energy 
eigenstates between the two glueball operators in C(t), it is easy to show 
that for large values of the euclidean time variable t, C(t) has the asymptotic 
behavior 

C{t) —>< vac\g{Q)\g >< g\g(0)\vac > exp(-m g t) + . . . , (6) 

where \g > is the zero-momentum state of the lightest scalar glueball, m g is its 
mass, and the omitted terms come from scalar glueball excitations and fall off 
exponentially in t with a coefficient larger than m g . So for large t, the effective 
glueball mass m g (t) has asymptotic behavior 

, \ , o(t) 

"»sW = log C (t + 1) ~* TOfi " ( ' 

Figure || shows m g (t) in lattice units, for a lattice 30 x 32 2 x 40 at f3 of 6.4, 
obtained from an ensemble of 25440 configurations of chromoelectric field. The 
lattice spacing in this case, determined from the p mass, is about 0.05 fm and 
lattice period is about 1.6 fm. For t of 4 and greater, m g (t) is consistent with 
a constant shown by the horizontal line, giving according to Eq. (^) a value for 
m g . From similar calculations for a range of different lattice periods, we found 
that 1.6 fm gives results within a fraction of a percent of the infinite- volume 
limit. 

To obtain the zero-lattice-spacing limit of m g , we evaluated C (t) for five 
different values the lattice spacing, in all cases with lattice period of 1.6 fm or 
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Figure 3: Continuum limit of the scalar glueball mass. 

greater. Now, the lattice version of the action Sa, defined in Eqs. (|l|), which 
occurs in a valence approximation calculation of C(t), replaces derivates by 
symmetric finite difference, obeying 

f(x + a)-f(x-a) _ df(x) | a 2 d 3 f(x) | 
2a dx 6 dx 3 

where a represents lattice spacing. The error in replacing an integral by a 
sum in Sa goes to zero more rapidly than a 2 . Thus for small a, Sa will be 
the continuum action plus an error of order a 2 . To leading order in a, m g 
measured in physical units will then also be its continuum value plus an error 
of order a 2 . Figure ^ shows the five values of m g in physical units as a function 
of a 2 in physical units. The unit A-^L in Figure || is the valence approximation 
A|jj, for which a value of 243.7 ±6.8 MeV waSj-pbtained in the course of 
our calculation of quark- antiquark meson masses u. To convert m g a and a 
to physical units, values of A-^jLet in lattice units were found by the two-loop 
Callan-Symanzik equation. 

As it turns out, the p mass in lattice units m p a scales almost perfectly with 
A^-La for the range of parameters in Figure |j| So Figure || can also be thought 
ofsimply as a plot of m g and a in units of m p but with the axes mislabeled 
by powers of 2A3MeV/m p . 

The three points at smallest lattice spacing in Figure || fit a straight line 
in a 2 quite well. Extrapolating to zero lattice spacing gives the limiting m g 
of 1740(71) MeV, in good agreement with the observed mass of /j(1710), 
as mentioned earlier in Section [j]. This result is rather insensitive to how 
the extrapolation is done. Even if we had arbitrarily taken either of the last 
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two points as the continuum limit, a procedure which is certainly less reliable 
than extrapolation, the answer would change by less than half of the 71 MeV 
statistical uncertainty. 



5 Decay Couplings 

We calculated coupling constants both for the decay of the lightest scalar 
glucball to pairs of pseudoscalars at rest and to pairs of pscudoscalars with 
oppositely directed momenta of magnitude 2ir/L, for lattice period L. For the 
sake of simplicity, I am only going to give details for the decay to pseudoscalars 
at rest. To measure the decay at rest we evaluated the three-point function 

C(t g ,t„) = <g(t g )J2^f(^)Tr f (0)>, (9) 
/ 

tt/0) = J d 3 xW u {x,t)j 5 X f ^ d (x,t), 

where g(t) is defined in Eqs. (g) and the A/ are an orthonormal set of SU(3) 
flavor matrices. The u, d and s quark masses were set equal, and chosen so 
that the energy of the lightest flavor-singlet, two-pseudoscalar, zero-momentum 
state coincides with the glueball mass m g . How to apply this calculation to 
the real world with a very different set of quark masses, I will explain in the 
next section. 

By inserting a complete set of energy eigenstates between g(t g ) and ir^{t v ) 
in C{t g ,t. n ) 1 you can show that at large t g — t w C{t g ,t^) has the asymptotic 
behavior 

C{t g ,U)^ (10) 

< vac\g(0)\irir >< -ktt\ ir^t^ir j (0)\vac > exp(—m g t g ) + . . . , 
/ 

where \nn > is the lightest zero- momentum, flavor-singlest state of two pseu- 
doscalars and the terms omitted fall off exponentially in t g with coefficients 
larger than m g . 

The coupling constant we are looking for can be found from < wac|g(0)|7r7r >. 



To obtain < vac|(7(0)|7r7r > from C(t g ,t n ) using Eq. (|10[), we need the value 
of < 7T7r| 2/ "/(^Tr)?!"/ (0)1 va > c >■ This factor carries the corrections to the 
decay calculation arising from 7r — 7r final state interactions. We extracted 
< 7T7r| 7Tj(t 7r )7rj(0)|wac > from the four-point function 

) =< 4(^) 7r /(°) 12 7r s(*T) 7r ff(°) > • ( u ) 

/ 9 
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The job of evaluating C(t„3, t^tni) and extracting < ttit\ J^f "/(O^r f(®)\ vac > 
I will not describe here. If we had simply taken for < irir\ wJt^Tr f(Q)\vac > 
its value assuming no tt — tt interaction, however, the final coupling constant 
would be changed by less than 15%. 

As a byproduct of finding the four-point function C(t V 3, ^2^1); we were 
also able to determine the first omitted term in Eq. (|l(j). To accelerate the 
approach to large t g — t n asymptopia, we subtracted this term from C^g,^), 
giving D(t g , ijr). The large t g — t w asymptotic behavior of D{t g , t^) is then 

D(tg,U) -► XK(tg,t„), (12) 

where A is the decay coupling constant for g — » rnr at rest and K(t g ,tn) 
is a kinematic factor determined in part from C(t n 3, t W 2, t^i)- Alternatively, 
defining from £)(t 9 ,f ff ) an effective X(t g ,t n ), 

M V'> = w$r r < 13 » 

it follows that for large t g ~t n , X(t g , t^) approaches the decay coupling constant 
A. 

Values of X(t g , t w ) were calculated on a lattice 16 3 x 24 at (J of 5.70, corre- 
sponding to lattice spacing about 0.15 fm and lattice period of about 2.3 fm, 
from an ensemble of 10500 configurations of A field. For t g — fixed, we found 
A(t ff ,t T ) independent of for > 3. Figure |i] shows X(t g ,t„), for t„ > 3, as 
a function of t g — t w . Largely as a consequence of the subtraction in defining 
D(t g ,i w ), the data is consistent with a constant from t g — t n of zero on out. 
This constant is the final value of A for glueball decay to pseudoscalars at rest. 



6 Comparison with Experiment 

So far I have given u, d and s quarks degenerate, unphysical mass values. 
Here is how to fix that. An expansion to first order in the quark mass matrix 
taken around some relatively heavy SU(3) symmetric point gives glueball decay 
couplings to 7r's, K's and the ry's which are a common linear function of each 
meson's average quark mass. Since meson masses squared are also nearly a lin- 
ear function of average quark mass, the decay couplings are a linear function 
of meson masses squared. Therefore, a linear fit to our predictions for de- 
cay couplings as a function of pseudoscalar mass squared lets you extrapolate 
from unphysical degenerate values of quark masses to physical nondegener- 
ate values of quark masses. Figure ^| shows predicted coupling constants as a 
function of predicted meson mass squared along with a linear extrapolation of 
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Figure 4: A as function of t g — t n . 



the predicted values to the physical n, K and 7/ masses. Shown also are the 
observed couplingfi for decays of /j(1710) to pairs of 7r's, K's and 77's. Every- 
thing is in units of the p mass. The total predicted width for glueball decay 
to pseudoscalar pairs becomes 108(28) MeV, in comparison to 99(15) MeV for 



How far are the valence approximation, finite lattice spacing decay cou- 
plings likely to be from the real world? From the comparison of finite lattice 
spacing valence approximation hadron masses with their values in the real 
world, I would expect an error of 15% or less in going to the continuum limit 
and another 6% or less arising from the valence approximation. The total 
predicted width for glueball decay to two pseudoscalars should then have an 
error of less than 50%. A 50% increase in our predicted two-body decay width, 
combined with any reasonable corresponding guess for multibody decays, gives 
a total glueball width small enough for the particle to be observed easily. 

For the continuum limit glueball mass, a 6% valence approximation error 
would be 100 MeV, but according to an adaptation oLan argument giving a 
negative sign for the valence approximation error in f^u, the sign of this error 
is also expected to be negative. Thus the the only established resonance aside 
from /j(1710) with the^iorrect quantum numbers and mass close enough to 
1740 to be a candidate □ for the scalar glueball is /o(1500). The most likely 
interpretation of /o(1500), however, I think is as an ss quark- antiquark meson. 
The su scalar and tensor are nearly degenerate at about 1430 MeV. So the ss 
scalar and tensor should lie close to each other somewhere above 1430 MeV. 
The ss tensor has been identified at 1525 MeV. An ss scalar around 1500 MeV 
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Figure 5: Decay couplings. 

seems to me hard to avoid. 
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